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We present the application of variational- wavelet analysis to numerical/analytical 
calculations of Wigner functions in (nonlinear) quasiclassical beam dynamics prob- 
lems. (Naive) deformation quantization and multiresolution representations are the 
key points. 



1 INTRODUCTION 

In this paper we consider some starting points in the applications of a new 
numerical-analytical technique which is based on local nonlinear harmonic 
analysis (wavelet analysis, generalized coherent states analysis) to the quan- 
tum/quasiclassical (nonlinear) beam/accelerator physics calculations. The 
reason for this treatment is that recently a number of problems appeared in 
which one needs take into account quantum properties of particles/beams. Our 
starting point is the general point of view of deformation quantization ap- 
proach at least on naive Moyal/Weyl/ Wigner level (part 2). The main point 
is that the algebras of quantum observables are the deformations of commuta- 
tive algebras of classical observables (functions) [l].So, if we have the Poisson 
manifold M (symplectic manifolds, Lie coalgebras, etc) as a model for clas- 
sical dynamics then for quantum calculations we need to find an associative 
(but non-commutative) star product * on the space of formal power series in 
h with coefficients in the space of smooth functions on M such that 



/ * g = fg + Hf, g} + Y, K n B n (f, g), 



(l) 



where {/, g} is the Poisson brackets, B n are bidifferential operators C°°(X) ® 
C°°(X) — > C°°(X). There is also an infinite-dimensional gauge group on the 
set of star-products 

(2) 



n>2 



where D n are differential operators. Kontsevich gave the solution to this de- 
formation problem in terms of formal power series via sum over graphs[l]. He 
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also proved that for every Poisson manifold M there is a canonically defined 
gauge equivalence class of star-products on M. Also there is the nonperturba- 
tive corrections to power series representation for * [1]. In naive calculations 
we may use simple formal rules: 

* = exp(^(*d x 'ff p -*dp^ x fj (3) 
f(x,p)*g(x,p) = f(x,p- yS,) -g(x,p + — d x ) (4) 
= f( x + y^p^P- l -7£~$x)g{x,p) (5) 

In this paper we consider calculations of Wigner functions (WF) as the 
solution of Wigner equations [2] (part 3): 

ih^-W(x,p,t) = H * f(x,p,t) - f(x,p,t) * H (6) 
and especially stationary Wigner equations: 



H*W = W*H = Ef (7) 

Our approach is based on extension of our variational-wavelet approach 
[3]- [14]. Wavelet analysis is some set of mathematical methods, which gives us 
the possibility to work with well-localized bases (Fig. 1) in functional spaces 
and gives maximum sparse forms for the general type of operators (differ- 
ential, integral, pscudodifferential) in such bases. These bases are natural 
generalization of standard coherent, squeezed, thermal squeezed states [2], 
which correspond to quadratical systems (pure linear dynamics) with Gaus- 
sian Wigner functions. So, we try to calculate quantum corrections to classical 
dynamics described by polynomial nonlinear Hamiltonians such as orbital mo- 
tion in storage rings, orbital dynamics in general multipolar fields etc. from 
papers [3]-[13]. The common point for classical/quantum calculations is that 
any solution which comes from full multiresolution expansion in all space/time 
(or phase space) scales represents expansion into a slow part and fast oscil- 
lating parts (part 4). So, we may move from the coarse scales of resolution to 
the finest one for obtaining more detailed information about our dynamical 
classical/quantum process. In this way we give contribution to our full solu- 
tion from each scale of resolution. The same is correct for the contribution to 
power spectral density (energy spectrum): we can take into account contribu- 
tions from each level/scale of resolution. Because affine group of translations 
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and dilations (or more general group, which acts on the space of solutions) is 
inside the approach (in wavelet case), this method resembles the action of a 
microscope. We have contribution to final result from each scale of resolution 
from the whole underlying infinite scale of spaces. In part 5 we consider nu- 
merical modelling of Wigncr functions which explicitly demonstrates quantum 
interference of generalized coherent states. 

2 Quasiclassical Evolution 

Let us consider classical and quantum dynamics in phase space f2 = R 2m 
with coordinates (x, £) and generated by Hamiltonian H(x, £) <G C°°(f2; R). If 
: £1 — ► Q is (classical) flow then time evolution of any bounded classical 
observable or symbol b(x,£) <G C°°(Sl,R) is given by bt(x,£) = b($^(a:,£)). 
Let H = Op w (H) and B — Op w (b) are the self-adjoint operators or quantum 
observables in L 2 (R n ), representing the Weyl quantization of the symbols H, b 
[1] 

{BU){X) = (2^ J R2n b ' * i<(x ~ vU>/ Mv)<iV<%, (8) 

where u £ S(R n ) and B t = e lHt / h Be~ iHt / n be the Heisenberg observable or 
quantum evolution of the observable B under unitary group generated by H. 
B t solves the Heisenberg equation of motion B t = (i/h)[H, B t \. Let bt(x, £; h) 
is a symbol of B t then we have the following equation for it 

bt = {H,b t } M , (9) 

with the initial condition bo(x, £ 5 K) — b{x, £). Here {/, q}m(x, £) is the Moyal 
brackets of the observables f,g& C°°(R 2n ), {/, g}tn{x, £) = f$g — g$f, where 
f$g is the symbol of the operator product and is presented by the composition 
of the symbols /, g 

(/«<?)(*, = j^f 2 j Rin e -<^>/^«^>/* (10) 
■f(x + w, p + £) ■ g(x + r, t + ^)dpdrdrdw 

For our problems it is useful that {/, <?}m admits the formal expansion in 
powers of K: 

{f,g}M{x,0~{f,g} + 2-i ■ (- 1 )' /31 • ( d U D *9) ■ (0fff££/),(ll) 

\a+0\=j>l 
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where a — (ai, . . . , a n ) is a multi-index, \a\ = a\ + . . . + a„, -Da; = —ihd x . 
So, evolution (1) for symbol b t (x,^; K) is 

6* = {W,6t} + ^ E (-l) 1 ^^^^).^^). (12) 

|a+/3|=j>l 

At H = this equation transforms to classical Liouville equation. Equation 
(12) plays the key role in many quantum (semiclassical) problems. We con- 
sider its particular case-Wigner equation-in the next section. 



3 Wigner Equations 

According to Weyl transform quantum state (wave function or density opera- 
tor) corresponds to Wigner function, which is analog of classical phase-space 
distribution [2]. We consider the following form of differential equations for 
time-dependent WF 



d t W(p, g, i) = - sin [- @f - d* ) ■ H(p, q)W(p, q, t) ( I A ) 



Let 



p = |* e >< * £ 



(14) 



be the density operator or projection operator corresponding to the en- 
ergy eigenstate |\l/ e > with energy eigenvalue e. Then time-independent 
Schroedinger equation corresponding to Hamiltonian 

(15) 



H(p,q) = ^- + U(q) 

where U(q is arbitrary polynomial function (related beam dynamics models 
considered in [3]-[13]) on q is [2]: 



Hp = ep 



(16) 



After Weyl- Wigner mapping we arrive at the following equation on WF in 
c-numbers: 



(17) 



or 



H p d Ti 2 d 2 



)w(p,q) + u(q-l-jp) W{p, q) = eW{p, q) 



(— + 

\2m 2imdq 8m dq 2 J \* 2idpy 
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After expanding the potential U into the Taylor series we have two real partial 
differential equations 

(_P_d_ ^ 1 (ihY m d 2m+1 U d 2m+1 \ 

\md~q + ^ (2m+l!) ItJ dq 2 ^ 1 Op 2 ^ 1 J ^ ?) ~ ( j 

m=0 v ' 

eW(p,?) (19) 

In the next section we consider variation-wavelet approach for the solution of 
these equations for the case of arbitrary polynomial U(q), which corresponds 
to a finite number of terms in equations (18), (19) up to any order of h. 

4 Variational Multiscale Representation 

Let L be arbitrary (non)linear differential operator with matrix di- 
mension d, which acts on some set of functions ^ = ^f(x,y) = 

(^>\x, y), * d (x, yfj, x,yenc$t 2 from L 2 (Q): 

L^ = L(Q,x,y)^(x,y) = 0, (20) 

where 

Q = Qd 1 42,d 3 ,dA x ,y,d/dx,d/dy) = ^ aijkiX%y3 idx) (~d~) ^ 



Let us consider now the N mode approximation for solution as the following 
ansatz (in the same way we may consider different ansatzes) : 



N 



* N (x,y)= a r , s * r (x)$s(y) (22) 

r,s — 1 

We shall determine coefficients of expansion from the following Galcrkin con- 
ditions (different related variational approaches are considered in [3]- [13]): 

C = j {L^ N )^ k {x)$> t {y)dxdy = (23) 

So, we have exactly dN 2 algebraical equations for dN 2 unknowns a rs . 

But in the case of equations for WF (18), (19) we have overdetermined 
system of equations: 27V 2 equations for N 2 unknowns a rs (in this case d = 1). 
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In this paper we consider non-standard method for resolving this problem, 
which is based on biorthogonal wavelet expansion. So, instead of expansion 
(22) we consider the following one: 




(24) 




where %?i(x)$j (y) are the bases dual to initial ones. Because wavelet functions 
are the generalization of coherent states we consider an expansion on this 
overcomplete set of bases wavelet functions as a generalization of standard 
coherent states expansion. 

So, variational/Galerkin approach reduced the initial problem (20) to the 
problem of solution of functional equations at the first stage and some alge- 
braical problems at the second stage. We'll consider now the multircsolution 
expansion as the second main part of our construction. Because affine group 
of translation and dilations is inside the approach, this method resembles the 
action of a microscope. We have contribution to final result from each scale 
of resolution from the whole infinite scale of increasing closed subspaces Vf. 
...V-2 C V-i C Vo C Vi C V2 C .... The solution is parametrized by solutions 
of two reduced algebraical problems, one is linear or nonlinear (23) (depends 
on the structure of operator L) and the second one are some linear problems 
related to computation of coefficients of algebraic equations (23). These co- 
efficients can be found by the method of Connection Coefficients (CC)[15] or 
related method [16]. We use compactly supported wavelet basis functions for 
expansions (22), (24). We may consider different types of wavelets including 
general wavelet packets (section 5 below). These coefficients depend on the 
wavelet-Galerkin integrals. In general we need to find (di > 0) 



According to CC method [15] we use the next construction for quadratic case. 
When N in scaling equation is a finite even positive integer the function ip(x) 
has compact support contained in [0,N — 1]. For a fixed triple (di,d,2,ds) 
only some A^f 2 ** 3 are nonzero: 2~N<i<N-2, 2-N<m<N- 
2, \£ - m\ < N - 2. There are M = 3N 2 - 9N + 7 such pairs (£,m). 
Let A dld2d3 be an M-vector, whose components are numbers A^ 2ds . Then 
we have the following reduced algebraical system : A satisfy the system of 
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00 




(25) 



— OO 



equations (d = di + d2 + d^) 

AA d ld2 d 3 = 2 i-d A d ld2 d^ At, m . q , r = a P a q -2i +P a r -2 m + P (26) 

p 

By moment equations we have created a system of M + d + 1 equations in 
M unknowns. It has rank M and we can obtain unique solution by combi- 
nation of LU decomposition and QR algorithm. For nonquadratic case we 
have analogously additional linear problems for objects (25). Solving these 
linear problems we obtain the coefficients of reduced main linear/nonlinear 
algebraical system (23) and after its solution we obtain the coefficients of 
wavelet expansion (22), (24). As a result we obtained the explicit solution of 
our problem in the base of compactly supported wavelets (22). 

Also in our case we need to consider the extension of this approach to 
the case of any type of variable coefficients (periodic, regular or singular). We 
can produce such approach if we add in our construction additional refinement 
equation, which encoded all information about variable coefficients [16]. So, 
we need to compute only additional integrals of the form 

b ij {t){i Pl ) d -(2 m t - fc 1 )(^) d2 (2 m t - fc 2 )dx, (27) 

D 

where bij(t) are arbitrary functions of time and trial functions ip\, ip% satisfy 
the refinement equations: 

W(t) = 5^0ifc^(2t-fe) (28) 
feez 

If we consider all computations in the class of compactly supported wavelets 
then only a finite number of coefficients do not vanish. To approximate the 
non-constant coefficients, we need choose a different refinable function 
along with some local approximation scheme 

{B t f){x) := Fe, k (f)M2 e t ~ fc), (29) 

where Fg t k are suitable functionals supported in a small neighborhood of 2~ e k 
and then replace bij in (27) by Bgbij(t). To guarantee sufficient accuracy of 
the resulting approximation to (27) it is important to have the flexibility of 
choosing ^3 different from (fi,(f2- So, if we take if4 = xd, where xd is 
characteristic function of D, which is again a refinable function, then the 



problem of computation of (27) is reduced to the problem of calculation of 
integral 



fcl) 
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ip 3 (2 e t - k 2 )<pt (2 r t - k 3 )^ 2 (2 s t - k 4 )dx 



(30) 



The key point is that these integrals also satisfy some sort of algebraical 
equation [16]: 

2-MH(k) = J2 b ^k-eH(£), M = <*i + d 2 . (31) 
iez 

This equation can be interpreted as the problem of computing an eigenvec- 
tor. Thus, the problem of extension of our approach to the case of variable 
coefficients is reduced to the same standard algebraical problem as in case of 
constant coefficients. So, the general scheme is the same one and we have only 
one more additional linear algebraic problem. After solution of these linear 
problems we can again compute coefficients of wavelet expansions (22), (24). 

Now we concentrate on the last additional problem which comes from 
overdeterminity of equations (18), (19), which demands to consider expansion 
(24) instead of expansion (22). It leads to equal number of equations and 
unknowns in reduced algebraical system of equations (23). For this reason 
we consider biorthogonal wavelet analysis. We started with two hierarchical 
sequences of approximations spaces [16]: . . . V_ 2 C V_i C V C V\ C V 2 ■ ■ ■ , 
. . . VI 2 C V_ 1 C Vo C V\ C V 2 . . . , and as usually, Wo is complement to Vo in 
V\, but now not necessarily orthogonal complement. New orthogonality con- 
ditions have now the following form: Wq _L Vo, Wq .1 Vo, Vj _L Wj, Vj _L 
Wj, translates of ip span Wo, translates of ip span Wo. Biorthogonal- 
ity conditions are < tpjk^j'k' >= J^° 00 ' l Pjk(x)ipj'k'{x)dx = 4fc'%', where 
ipjk{x) = 2 J / 2 i/)(2 J 'x — k). Functions tp(x),ip(x — k) form dual pair: < tp(x — 
k),ip(x — £)>=5ki, < <f{x — k),ip(x — £) >— for Vfc, \/£. Functions ip, (p 
generate a multiresolution analysis. ip(x — k), ip(x — k) are synthesis functions, 
ip(x — £), xp(x — £) are analysis functions. Synthesis functions are biorthogonal 
to analysis functions. Scaling spaces are orthogonal to dual wavelet spaces. 
Two multiresolutions are intertwining Vj + Wj = Vj+i, Vj + Wj — Vj+i. 
These are direct sums but not orthogonal sums. So, our representation for 
solution has now the form f(t) = J2j k bjkipjk{t), where synthesis wavelets are 
used to synthesize the function. But bjk come from inner products with anal- 
ysis wavelets. Biorthogonal point of view is more flexible and stable under the 
action of large class of operators while orthogonal (one scale for multiresolu- 
tion) is fragile, all computations are much more simple and we accelerate the 
rate of convergence of our expansions (24). By analogous anzatzes and ap- 
proaches we may construct also the multiscale/multiresolution representations 
for solution of time dependent Wigner equation (13) [14]. 
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Figure 1. Localized contributions to beam motion. 



5 Numerical Modelling 

So, our constructions give us the following N-mode representation for solution 
of Wigner equations (18)-(19): 

N 

W N (p,q)= a rs * r (p)$s(q) (32) 

r,s — 1 

where ^ r {p), ®s(q) may be represented by some family of (nonlinear) eigen- 
modes with the corresponding multiresolution/multiscale representation in 
the high-localized wavelet bases (Fig. 1): 

i>Mi 

i>M 2 

Our (nonlinear) eigenmodes are more realistic for the modelling of nonlinear 
classical/quantum dynamical process than the corresponding linear gaussian- 
like coherent states. Here we mention only the best convergence properties 
of expansions based on wavelet packets, which realize the so called minimal 
Shannon entropy property (Fig. 1). On Fig. 2 we present numerical modelling 
[17] of Wigner function for a simple model of beam motion, which explicitly 
demonstrates quantum interference property. On Fig. 3 we present the multi- 
scale/multiresolution representation (32)-(34) for solution of Wigner equation. 
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Signal in time 




WV, lin. scale. Tnreshold=5% 




Frequency [Hz] Timfj [s] 

Figure 2. Wigner function for 3 wavelet packets. 
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Figure 3. Multircsolution/multiscale representations for Wigner functions. 
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